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Abstract 



It is shown that a top mass of 174 ± 17 GeV, as quoted recently by the CDF Col- 
laboration, constrains the mixing angle between the Higgs doublets in the Minimal 
Supersymmmetric extension of the Standard Model (MSSM) to: 1.2 < tan/3 < 5.5 
at the 90% C.L.. The most probable value corresponds to tan/3 = 1.56; such a small 
value causes a large mixing in the stop sector and the lightest stop is likely to be 
below the top mass. In this case the stop production in pp collisions would contribute 
to the top signature, thus providing a possible explanation for the large effective ti 
cross section observed by CDF. 
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1 Introduction 



The Grand Unification idea has been subjected recently to a new test using the new precise 
LEP data [|l], |2], M. The result clearly indicates that the minimal Standard Model (SM) 
does not lead to unification of the coupling constants, if they are extrapolated to high en- 
ergies |2j. On the contrary, within the Minimal Supersymmetric extension of the Standard 
Model (MSSM) unification is achieved. Supersymmetry [|j presupposes a symmetry between 
fermions and bosons, thus doubling the particle spectrum of the SM. The predicted particles 
are indicated by a tilde above the usual SM symbol. Since these supersymmetric particles 
( "sparticles" ) have not yet been observed, supersymmetry must be broken. From the unifi- 
cation condition a first estimate of the SUSY breaking scale could be made: it was found to 
be of the order of 1000 GeV, or more precisely 10 3±1 GeV 0]. 

Assuming soft symmetry breaking at the Grand Unification (GUT) scale, all sparticle 
masses can be expressed in terms of 5 parameters and the masses at low energies are then 
determined by the well known Renormalization Group (RG) equations. The parameters are: 
mo, the common mass of the spin squarks and sleptons; my 2 , the common mass of the 
spin 1/2 gauginos; /i, the mixing parameter between the Higgs doublets; tan/3, the ratio of 
the vacuum expectation values of the two Higgs doublets; and A, the trilinear coupling in the 
Higgs sector. So many parameters cannot be derived from the unification condition alone. 
Further constraints can be considered: 



Mz predicted from electroweak symmetry breaking p] - [13]. 



Constraints from the unification of Yukawa couplings [10, 11], [14] 
Constraints from the lower limit on the proton lifetime [24, ^5|, 26] 
Experimental lower limits on SUSY masses [27, 28]. 



• Constraints from the top mass suggested by CDF [29|. 

We perform a statistical analysis, in which all constraints are implemented in a x 2 defi- 
nition and try to find the most probable region of the parameter space by minimizing the x 2 
function. The results will be presented after a short description of the experimental input 
values and it is shown that a likely solution has the lightest stop mass below the top mass. 



2 Unification of the Couplings 

In the SM based on the group SU(3) x SU(2) x U(l) the couplings are defined as: 

a x = (5/3)£/ ,2 /(4tt) = 5a/(3cos 2 ^) 

a 2 = # 2 /(47r) = a/sm 2 6 w (1) 

as = 5s/(4vr) 

where g' ,g and g s are the U(l), SU(2) and SU(3) coupling constants; the first two coupling 
constants are related to the fine structure constant by: e = V4vra = gsinOw = g' cos 9w- 

In the MS renormalization scheme the world averaged values of the coupling constants 
at the Z° energy are 

aT 1 {M z ) = 127.9 ±0.1 (2) 
sin 2 ^ = 0.2324 ± 0.0005 (3) 
a 3 = 0.123 ±0.006. (4) 



The value of a~ l is given in ref. [^] and the value of sin 2 Ojfg has been been taken from a 
detailed analysis of all available data by Langacker and Polonsky ]3l| , which agrees with the 
latest analysis of the LEP data ||32|| . The error includes the uncertainty from the top quark. 
We have not used the smaller error of 0.0003 for a given value of mt, since the fit was only 
done within the SM, not the MSSM, so we prefer to use the more conservative error including 
the uncertainty from mt- 

The «3 value corresponds to the value at Mz as determined from quantities calculated 



in the "Next to Leading Log Approximation" [ 33 1 . These quantities are less sensitive to the 



renormalization scale, which is an indicator of the unknown higher order corrections; they 
are the dominant uncertainties in quantities relying on second order QCD calculations [34]. 



This a s value is in excellent agreement with a preliminary value of 0.120 ±0.006 from a fit to 
the Z° cross sections and asymmetries measured at LEP |32[| , for which the third order QCD 
corrections have been calculated too; the renormalization scale uncertainty is correspondingly 
small. 

The top quark mass was simultaneously fitted to all electroweak data and found to be Q| : 

M top = im±H -11 GeV > (5) 

where the first error is statistical and the second error corresponds to a variation of the 
Higgs mass between 60 and 1000 GeV. The central value corresponds to a Higgs mass of 300 
GeV. Preliminary analysis including the LEP 1993 data find M top = 165 ± 12 ± 18 GeV and 
M top = 174 ± llt\l GeV > if the 

new SLD data from SLAC is included |35(| . These values are 



in good agreement with recent results quoted by the CDF Collaboration J29[ |: 

M top = 174ijg t\l GeV, (6) 

where the first error is statistical and the second error systematic. 

For SUSY models, the dimensional reduction DR scheme is a more appropriate renormal- 
ization scheme p6[| . This scheme also has the advantage that all thresholds can be treated 
by simple step approximations. Thus unification occurs in the DR scheme if all three a^ 1 (fi) 
meet exactly at one point. This crossing point then gives the mass of the heavy gauge bosons. 
The MS and DR couplings differ by a small offset 

i i a 



(7) 



a DR a MS Un 

where the C, are the quadratic Casimir coefficients of the group (Cj = N for SU(iV) and 
for U(l) so ai stays the same). Throughout the following, we use the DR scheme for the 
MSSM. 



3 Mz from Electroweak Symmetry Breaking 

In the MSSM at least two Higgs doublets have to be introduced. Radiative corrections from 
the heavy top and stop quarks can drive one of the Higgs masses negative, thus causing 
spontaneous symmetry breaking in the electroweak sector. In this case the Higgs potential 
does not have its minimum for all fields equal zero, but the minimum is obtained for non- 
zero vacuum expectation values of the fields. The scale, where symmetry breaking occurs 
depends on the starting values of the mass parameters at the GUT scale, the top mass and 
the evolution of the couplings and masses. This gives strong constraints between the known 
mass and the SUSY mass parameters, as demonstrated e.g. in ref. JlOj ] . 



After including the one- loop corrections to the potential [37, BM, |3^, EO, 41, 20 1, the Mz 



mass becomes dependent on the top- and stop quark masses too. The corrections are zero if 
the top- and stop quark masses are identical, i.e. if supersymmetry would be exact. They 
grow with the difference mf — mf, so these corrections become unnaturally large for large 



values of the stop masses [42 



4 mi, from the m&/m r Mass Ratio 

Unification of the Yukawa couplings for a given generation at the GUT scale predicts relations 
for quark and lepton masses within a given family. Unfortunately, for the light quarks the 
masses are uncertain, but the ratio of b-quark and r-lepton masses can be correctly predicted 
by the radiative mass corrections @, |, ||, 0]. 

Assuming the simplest possible GUT model based on SU(5) gauge group, one has at the 
GUT scale: rrib = rn T . To calculate the experimentally observed mass ratio the RG equations 
for the running masses have to be used. By a physical mass we understand the value of the 
running mass at the energy scale equal to the mass itself. This definition of the mass is used 
throughout this paper. 

From the RG equations for the Yukawa couplings one can easily obtain the RGE for the 



ratio |42 




R b r^ = 

For the running mass of the b-quark we used 

m b = 4.25 ± 0.3 GeV. (8) 
This mass depends on the choice of scale and the value of a 3 {m^). Consequently, we have 



assigned a rather conservative error of 0.3 GeV instead of the proposed value of 0.1 GeV [45|. 



pole 



Note that the running mass (in the MS scheme) is related to the physical (pole) mass 

m b = Mf e 12.4(^) 2 ) « 0.825 M p b ° le , (9) 



by @ : 



so mj, = 4.25 corresponds to M b ° le ~ 5 GeV. We ignore the running of m T below mj, and use 



for the r mass: M T = 1.7771 ± 0.0005 GeV |6[. 



5 Top Mass Constraints 

The top mass can be expressed as: 

m t 2 = (4vr) 2 Y t {t) v 2 sin 2 (/3), (10) 



where the running of the Yukawa coupling Yt as function of t = log( ^y T ) in first ordery is 
given by @ : 

Y t (0)E(t) 
Yt(t) ~ l + 6y t (0)F(t)' 



4 Throughout the analysis we have used the second order RG equations, for which no analytical solution 
exists, but this will not change the following arguments dramatically. 



where E and F are functions of the couplings only (see appendix). One observes that Yt(t) 
becomes independent of Yt(0) for large values of Yt(0), implying an upper limit on the top 
mass. Requiring electroweak symmetry breaking implies a minimal value of the top Yukawa 
coupling, typically Y t (0) > 0(1CT 2 ). In this case the term 6Yt(0)F(t) in the denominator of 
eq. ( Jll| ) is much larger than one, since F(t) ~ 290 at the weak scale, where t ~ 66. In this 
case Y t (t) = E(t)/6F(t), so from eq. © it follows: 

ml = v 2 sin 2 (/3) » (190 GeV) 2 sin 2 (/3), (12) 

The physical (pole) mass is about 6% larger than the running mass : 

Mf oIe = m t (l + « (200 GeV) sin/3, . (13) 

The electroweak breaking conditions require 7r/4</3<7r/2; hence the equation above 
implies for the MSSM approximately: 

145 < Mf ole < 200 GeV, (14) 

which is consistent with the experimental values given in eqns. (|) and (|). 

For large top masses, the b-quark mass becomes a sensitive function of nit and of the 



starting values of the gauge couplings at Mqut [11 1- 



6 Experimental Lower Limits on SUSY Masses 

SUSY particles have not been found so far and from the searches at LEP one knows that the 
lower limit on the charged leptons and charginos is about half the Z° mass (45 GeV) p8[ 
and the Higgs mass has to be above 62 GeV (27|. The lower limit on the lightest neutralino 
is 18.4 GeV [28], while the sneutrinos have to be above 41 GeV [28]. These limits require 
minimal values for the SUSY mass parameters. 

There exist also limits on squark and gluino masses from the hadron colliders [28], but 
these limits depend on the assumed decay modes. Furthermore, if one takes the limits given 
above into account, the constraints from the limits of all other particles are usually fulfilled, 
so they do not provide additional reductions of the parameter space in case of the minimal 
SUSY model. 



7 Proton Lifetime Limits 

GUT's predict proton decay and the present lower limits on the proton lifetime yield quite 
strong constraints on the GUT scale and the SUSY parameters. The direct decay p — ► e + ir° 
via s-channel exchange requires the GUT scale to be above 10 15 GeV. This is not fulfilled 
in the SM, but always fulfilled in the MSSM. Therefore we do not consider this constraint. 
However, the decay via box diagrams with winos and Higgsinos predict much shorter lifetimes, 
especially in the preferred mode p — ► VK + . From the present experimental lower limit of 
10 32 yr for this decay mode Arnowitt and Nath [25] deduce an upper limit on the parameter 
B: 

B < (293 ± 42) M H J3M GVT GeV- 1 (15) 



Here Mu z is the Higgs triplet mass, which is expected to be of the order of Mqut- To obtain 
a conservative upper limit on B, we allow Mh 3 to become an order of magnitude heavier 
than Mqut, so we require 

B < 977 ± 140 GeV^ 1 . (16) 

The uncertainties from the unknown heavy Higgs mass are large compared with the contri- 
butions from the first and third generation, which contribute through the mixing in the CKM 
matrix. Therefore we only consider the second order generation contribution, which can be 



written as [25]: 



B = ^hs^™ 6 (17) 
as sm(2p) m~ 

where ct2 and 03 are the coupling constant of the SU(2) and SU(3) groups at the SUSY 
scale, respectively. One observes that the upper limit on B favours small gluino masses nig, 
large squark masses rrig, and small values of tan j3. To fulfill this constraint requires 

tan/? < 10 (18) 



for practically the whole parameter space [25 



8 Fit Strategy 

Prom the five parameters in the MSSM plus the common coupling acuT at the unification 
scale -Mgut one can determine all other SUSY masses, the b-quark mass, and Mz by per- 
forming the complete evolution of the couplings and masses including all thresholds. Details 
can be found in gg]. 

The most probable parameter values were obtained by minimizing the following x 2 func- 
tional: 
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5 We use the MINUIT program from F. James and M. Roos, MINUIT Function Minimization and Error 
Analysis, CERN Program Library Long Writeup D506; Release 92.1, from March 1992. Our \ 2 nas disconti- 
nuities due to the experimental bounds on various quanitities, which become "active" only for specific regions 
of the parameter space. Consequently the derivatives are not everywhere defined. The option SIMPLEX, 
which does not rely on derivatives, can be used to find the monotonous region and the option MIGRAD to 
optimize inside this region. 



The first term is the contribution of the difference between the three calculated and measured 
gauge coupling constants at Mz and the following two terms are the contributions from the 
Mz -mass and m\, -mass constraints. The last three terms impose constraints from the proton 
lifetime limits, from electroweak symmetry breaking, i.e. D = Vh(vi,V2) — Vff(0, 0) < 0, and 
from experimental lower limits on the SUSY masses. The top mass, or equivalently, the top 
Yukawa coupling enters sensitively into the calculation of mi, and Mz- Instead of the top 
Yukawa coupling one could have taken the top mass as a parameter. However, if the couplings 
are evolved from Mqijt downwards, it is more convenient to run also the Yukawa coupling 
downward, since the RG equations of the gauge and Yukawa couplings form a set of coupled 
differential equations in second order. Once the Yukawa coupling is known at Mgut, the top 
mass can be calculated at any scale. 

The following errors were attributed: cTj are the experimental errors in the coupling 
constants, as given above, a"b=0.3 GeV, <jb = 140 GeV -1 , while <jd and were set to 10 
GeV. The values of the latter errors are not critical, since the corresponding terms in the 
numerator are zero in case of a good fit and even for the 90% C.L. limits these constraints 
could be fulfilled and the x 2 was determined by the other terms, for which one knows the 
errors. 

The light thresholds are taken into account in the evolution of the coupling constants by 
changing the coefficients of the RGE at the value Q = rrii, where the threshold masses m« are 
obtained from the analytical solutions of the corresponding RGE. These solutions depend on 
the integration range, which was chosen between m; and Mgut- However, since one does 
not know mi at the beginning, an iterative procedure has to be used: one first uses Mz as a 
lower integration limit, calculates mj, and uses this as lower limit in the next iteration. Of 
course, since the coupling constants are running, the latter have to be iterated too, so the 
values of an(mi) have to be used for calculating the mass at the scale mj f|7]]. Usually 
three to five iterations are enough to find a stable solution. 

Following Ellis, Kelley and Nanopoulos [48| the possible effects from heavy thresholds are 
set to zero, since the proton lifetime limits forbid the Higgs triplet masses to be below Mgut • 
These heavy thresholds have been considered by other authors for different assumptions [49, 
0,0- 



9 Results 

We first consider fits without proton lifetime constraints, since they are only important in the 
determination of lower limits, as will be discussed below. The upper part of fig. [l] shows the 
evolution of the coupling constants in the MSSM for two cases: one for the minimum value 
of the x 2 function given in eq. |l9| (solid lines) and one corresponding to the 90% C.L. upper 
limit of the thresholds of the light SUSY particles (dashed lines). The position of the light 
thresholds is shown in the bottom part as jumps in the first order (3 coefficients, which are 
increased as soon as a new threshold is passed. Also the second order coefficients are changed 
correspondingly, but their effect on the evolution is not visible in the top figure in contrast to 
the first order effects, which change the slope of the lines considerably in the top figure. One 
observes that the changes in the coupling constants occur in a rather narrow energy regime, 
so qualitatively this picture is very similar to the case, in which all sparticles were assumed to 
be degenerate at an effective SUSY mass scale Mstjsy @- Since the running of the couplings 
depends only logarithmically on the sparticle masses, the 90% C.L. upper limits are as large 
as several TeV, as shown by the dashed lines in fig. [l] and more quantitatively in table |l|. 



With the fitted SUSY parameters given at the top of the table, the corresponding masses of 
the SUSY particles can be calculated. Their values are given in the lower part of the table. 

The 90% C.L. upper and lower limits on the masses are obtained by scanning mo and 
mi/2 till the x 2 value increases by 1.64, while optimizing the values of tan /3, //, , acuT; ^t(O) 
and Mqut- The lower limits on the SUSY parameters are shown in the left column of table |l|. 
The lowest values of mo and m^2 are required to have simultaneously a sneutrino mass above 
42 GeV and a wino mass above 45 GeV. If the proton lifetime limit is included, either mo or 
mi j2 have to be above a certain limit (see eq. [17]). Since the squarks and gauginos are much 
more sensitive to m.1/2 than mo, one obtains the lower limits by increasing mo- The minimum 
value for mo is about 400 GeV in this case. But in both cases the x 2 increase for the lower 
limits is due to the b-mass, which is predicted to be 4.6 GeV from the parameters determining 
the lower limits, so it gives a contributions to the x 2 function, which requires m;,=4.25 GeV 
(see eq. |8|). The 90% C.L. upper limits can be several TeV. If one requires that only solutions 
are allowed for which the corrections to Mz are not large compared to Mz itself, one has 
to limit the mass of the heaviest stop quark to about one TeV. The corresponding 90% C.L. 
upper limits of the individual sparticles masses are given in the right hand column of table 

The correction to Mz is 6 times Mz in this case. 

The mass of the lightest Higgs particle, called h in table [ll is a rather strong function 
of nit i as shown in fig. Q for various choices of tan/3, mo and m 1 / 2 - All other parameters 
were optimized for these inputs and after the fit the values of the Higgs and top mass were 
calculated and plotted. One observes that the mass of the lightest Higgs particle varies 
between 60 and 150 GeV and the top mass between 134 and 190 GeV. Furthermore, it is 



evident that tan/3 almost uniquely determines the value of mt (through eq. 1C), since even 
if nil 12 an d mo are varied between 100 and 1000 GeV, one finds practically the same m t for 
a given tan/3. The value of mt varies between 134 and 190 GeV, if tan/3 is varied between 



1.2 and 10. This range is in excellent agreement with the estimates given in eq. |Lj, if one 
takes into account that Mf ole ~ 1.06mj (see eq. 13). The shaded area in fig. ^ indicates the 



results on the top mass quoted by the CDF Collaboration [29]. It clearly favours low values 
of tan/3. Adding to the x 2 a term (Mt — 174) 2 /17 2 yields after minimization: 

1.2 < tan (3 < 5.5 at the 90% C.L. (20) 

The most probable value corresponds to tan/3 = 1.56 as indicated by the star in the figure. 

Such a low value leads to a large mixing in the stop sector, in which case a likely value of 
the lightest stop is below the top mass (see the typical fit in table [l]) , although stop masses 
above the top mass are not excluded, as shown by the upper limits in table ffl. Also a change 
in sign of the Higgs mixing parameter fi leads to stop masses above the top mass, but the x 2 
value is hardly worse in that case, so this cannot be excluded either. Varying -At(O) between 
+3mo and — 3mo does not influence the results very much, since it is usually compensated 
by a change in fx, so At(0) was kept zero, but its non-zero value at lower energies was taken 
into account. 

If the stop mass is below the top mass, it cannot decay into the top, but can decay as 
follows: 

h xt + b -» Xi + W + b -> x? + lepton + v + b, 



which is experimentally very similar to the normal top decay signature [51]. Additional stop 
production could be an explanation for the excess of events seen by the CDF Collaboration: 
they observe an effective cross section for top pair production of 13. 91^ pb, while the 
calculated tt cross section is only 5.8^ ' 4 pb [29]. 
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Symbol 


Lower limits 


Typical fit 


90% C.L. Upper limits 


Constraints 


GEY 


GEY+P 


GEY+(PF) 


GEY+ (P) 


GEY+(P)+F 


Fitted SUSY parameters 


m 


65 


400 


400 


400 


400 


m l/2 

1 


37 


80 


111 


1600 


475 




-117 


330 


870 


1842 


1101 


tan/3 


3.0 


3.0 


1.56 


8.5 


2.9 


Yt(0) 


0.0158 


0.0035 


0.0150 


0.0023 


0.0084 


M pole 






175 


178 


189 


m t 






165 


168 


178 


1/ogut 


23.8 


24.3 


24.5 


25.9 


25.2 


M G ut 


2.3 10 16 


2.0 10 16 


2.0 10 16 


0.8 10 16 


1.3 10 16 


SUSY masses in [GeV] 


X?(7) 


18 


25 


41 


720 


202 


X 2 (Z) 


39 


52 


80 


1346 


386 


xf(w) 


46 


48 


79 


1347 


386 


9 


109 


217 


293 


3377 


1105 


eL 


82 


406 


409 


1160 


521 


en. 


67 


401 


402 


729 


440 


VL 


41 


400 


406 


1157 


516 




120 


443 


477 


3030 


1071 




115 


440 


471 


2872 


1030 


~b L 


112 


352 


369 


2610 


903 


b R 


119 


440 


471 


2862 


1027 


h 






144 


2333 


725 








467 


2817 


1008 


x° 3 m 


109 


292 


540 


1771 


799 


xim 


120 


313 


556 


1780 


812 


xt(H±) 


129 


315 


566 


1816 


831 


h 






87 


146 


127 


H 


118 


523 


812 


2218 


1033 


A 


92 


521 


810 


2217 


1031 


H± 


121 


527 


813 


2219 


1034 



Table 1: Values of SUSY masses and parameters for various constraints: G=gauge coupling 
unification; E=electroweak symmetry breaking; Y= Yukawa coupling unification; P=Proton 
lifetime constraint; F=finetuning constraint. Constraints in brackets indicate that they are 
fulfilled but not required. The minimum values of the lightest Higgs mass, the stop mass and 
the top mass can't be reached for the parameters minimizing the squarks and slepton masses. 
One needs smaller values of tan (3 in that case. 
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Figure 1: Evolution of the inverse of the three couplings in the MSSM. The line above 
Mqut follows the prediction from the supersymmetric SU(5) model. The SUSY thresholds 
have been indicated in the lower part of the curve: they are treated as step functions in the 
first order (3 coefficients in the renormalization group equations, which correspond to a change 
in slope in the evolution of the couplings in the top figure. The dashed lines correspond to 
the 90% C.L. upper limit for the SUSY thresholds. 




Figure 2: The mass of the lightest Higgs particle as function of the top quark mass for values 
of tan/3 between 1.2 and 10 and values of mo and m\i2 between 100 and 1000 GeV. The 
parameters of fj,, Mqtjt, «gut and V*(0) are optimized for each choice of these parameters; 
the corresponding values of the top and lightest Higgs mass are shown as symbols. For 
small values of the Higgs mass increases with mo, as shown for a "string" of points, 
each representing a step of 100 GeV in mo for a given value of m 1( / 2 , which is increasing 
in steps of 100 GeV, starting with the low values for the lowest strings. At high values 
of mi/2 the value of mo becomes irrelevant and the "string" shrinks to a point. Note the 
strong positive correlation between mhiggs and all other parameters: the highest value of the 
Higgs mass corresponds to the maximum values of the input parameters, i.e. tan (5 = 10, 
mo = mi/2 = 1000 GeV; this value does not correspond to the minimum x 2 . More likely 
values are: m higgs r* 87 GeV for m 1/2 = 100 GeV, m = 400 GeV, \i = 822 GeV and 
tan/3 = 1.6, as indicated by the star. The hatched area corresponds to the top mass range 
measured by [29]. 
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